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X-RAY EMISSION FROM THE PLANET PULSAR B1257+12 
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ABSTRACT 

We report the detection of the millisecond pulsar B1257+12 with the Chandra X-ray Observatory. 
In a 20 ks exposure we detected 25 photons from the pulsar, with energies between 0.4 and 2.0 keV, 
corresponding to the flux Fx — (4.4 ± 0.9) x 10~ 15 ergs s" 1 cm~ 2 in this energy range. The X- 
ray spectrum can be described by a power-law model with photon index T « 2.8 and luminosity 
Lx ~ 2.5 x 10 29 ergs s _1 in the 0.3-8 keV band, for a plausible distance of 500 pc and hydrogen 
column density Ah = 3 x 10 20 cm -2 . Alternatively, the spectrum can be fitted by a blackbody model 
with kT rs 0.22 keV and projected emitting area ~ 2000 m 2 . If the thermal X-rays are emitted from 
two symmetric polar caps, the bolometric luminosity of the two caps is 2L D oi ~ 3 x 10 29 ergs s _1 . 
We compared our results with the data on other 30 millisecond pulsars observed in X-rays and found 
that the apparent X-ray efficiency of PSR B1257+12, Lx/E ~3x 10~ 5 for d = 500 pc, is lower than 
those of most of millisecond pulsars. This might be explained by an unfavorable orientation of the 
X-ray pulsar beam if the radiation is magnetospheric, or by strong asymmetry of polar caps if the 
radiation is thermal (e.g., one of the polar caps is much brighter than the other and remains invisible 
for most part of the pulsar period). Alternatively, it could be attributed to absorption of X-rays in 
circumpulsar matter, such as a flaring debris disk left over after formation of the planetary system 
around the pulsar. 

Subject headings: pulsars: individual (PSR B1257+12 = J1300+1240, PSR B1620-26 = J1623-2631) 
— stars: neutron — planetary systems — X-rays: stars 



1. INTRODUCTION 

Most of millisecond radio pulsars (MSPs) reside in bi- 
nary systems, usually with a white dwarf companion. 
It is generally believed that they have been spun up 
to millisecond periods by accretion from their compan- 
ions. About 20% of the known MSPs are solitary pulsars 
(Lorimer 2005 1 ), which either lost their binary compan- 
ions or were born with short periods and low magnetic 
fields. Both binary and solitary millisecond pulsars can 
emit X-rays from their magnetospheres. Similar to "ordi- 
nary" (non-recycled) pulsars, this magnetospheric X-ray 
emission is usually characterized by a power-law (PL) 
spectrum with photon indices T w 1-2 and sharp X-ray 
pulsations (e.g., Becker & Pavlov 2001; Zavlin 2006ab). 
In addition to the magnetospheric emission, MSPs emit 
thermal X-rays from hot polar caps. The polar cap emis- 
sion is characterized by a soft X-ray spectrum, which re- 
sembles a blackbody (BB) spectrum with temperatures 
~ 1-3 MK, and smoother pulsations. Both thermal and 
magnetospheric luminosities are small fractions, ~ 10~ 4 - 
10~ 2 ' 5 , of the spin-down power E; the thermal emission 
apparently dominates in solitary MSPs with E < 10 35 
ergs s _1 . 

Of particular interest among the solitary millisec- 
ond pulsars is PSR B1257+12 (also known as PSR 
J1300+1240; we will call it B1257 hereafter), which hosts 
a planetary system comprised of the first extrasolar plan- 
ets discovered (Wolszczan & Frail 1992). The pulsar's 
period, period derivative, and dispersion measure are 
P = 6.22 ms, P obs = 1.14 x 10~ 19 s s"\ and DM = 10.2 
cm -3 pc, respectively. The DM-based distance estimates 

Electronic address: pavlov@astro.psu.edu 

1 See http://www.livingreviews.org/lrr-2005-7 



are 620 and 450 pc for the Galactic electron density dis- 
tibution models by Taylor & Cordes (1993) and Cordes 
& Lazio (2002), respectively, with nominal uncertainties 
of 20%. Wolszczan et al. (2000) reported a distance of 
~800 pc from a timing parallax measurement, but this 
value is subject to large errors because of strong corre- 
lations with the other timing parameters. We will scale 
the distance to 500 pc below. 

B1257 shows a substantial proper motion, 46.4 ± 0.1 
and —82.2 ± 0.2 mas yr _1 in the right ascension and 
declination, respectively (Wolszczan et al. 2000), corre- 
sponding to the transverse velocity, v± = 224 d 500 km 
s , unusually large for a millisecond pulsar. There- 
fore, the period derivative should be corrected for the 
effect of pulsar's transverse motion (Shklovskii 1970): 

P = -Pobs - {vlP/dc) = (11.42 - 6.74 d 5 oo) x lO" 20 . 
This distance-dependent correction increases the inferred 
spin-down age, r = 0.863 (1 — 0.590d 50 o) _1 Gyr, and re- 
duces the inferred spin-down power, E = 1.88 x 10 34 (1 — 
0.590c?5oo) ergs s -1 , and magnetic field, B = 8.53 x 
10 8 (1 - 0.590d 50 o) 1/2 G. Notice that allowance for this 
effect puts an upper limit d < 848 pc, on the the distance 
to B1257. 

One can expect that B1257, just as other MSPs, emits 
X-rays generated in the pulsar's magnetosphere and/or 
hot polar caps, and studying its X-ray emission can be 
used, together with the data on other MSPs, to under- 
stand how the polar cap and magnetosphere properties 
depend on pulsar parameters. On the other hand, it is 
possible that some kind of debris in the pulsar's plane- 
tary system can lead to additional absorption of X-rays 
close to the pulsar, which might allow one to establish 
the presence of such debris and study their properties. 




Fig. 1.— ACIS-S3 image of B1257+12. The cross marks the radio 
position of PSR B1257+12, whose uncertainty is negligible on this 
scale. The 0'.'6 radius of the smaller (blue) circle corresponds to 
the 90% uncertainty of the Chandra pointing. The larger (white) 
circle of 1'.'5 radius encompasses the area from which the photons 
used for the spectral analysis were extracted. 

With this in mind, we proposed to observe B1257 with 
the Chandra X-ray Observatory. We describe the ob- 
servation and its results in §2 and §3, and discuss some 
implications in §4. 

2. OBSERVATIONS 

We observed B1257 with the Advanced CCD Imaging 
Spectrometer (ACIS) aboard Chandra on 2005 May 22 
(start time 53,512.220183 MJD) for 20.05 ks (19,797 s 
effective exposure time). The observation was carried 
out in Very Faint mode, and the pulsar was imaged on 
ACIS-S3 chip with a standard Y offset of -0'.'33. The de- 
tector was operated in Full Frame mode, which provides 
time resolution of 3.2 seconds. The data were reduced 
using the Chandra Interactive Analysis of Observations 
(CIAO) software (ver. 3.3; CALDB 3.2.1). For the anal- 
ysis, we used the standard grade filtering and restricted 
the energy range to 0.3-8 keV. 

3. X-RAY IMAGE AND SPECTRUM 

Figure 1 shows the ACIS-S3 image of the field around 
B1257, with an X-ray source centered at R.A. = 
13 h 00 m 03?10, deck = +12°40'56"0 (J2000). The uncer- 
tainty of this position, 0'.'6 at the 90% confidence level, 
is mainly caused by errors in absolute Chandra astrom- 
etry. Since it differs by only 0'.'3 from the radio position, 
R.A. = 13 h 00 m 03?0810, deck = +12°40'55'/875 for the 
epoch of the Chandra observation, we conclude with con- 
fidence that we detected the X-ray emission from B1257. 
The distribution of source counts in the ACIS image is 
consistent with that of a point source. 

To measure the flux and the spectrum of the pulsar, we 
chose a circular aperture of l'/5 radius (about 3 ACIS pix- 
els), which contains 95% encircled energy fraction. Using 
the CIAO psextract task, we found 25 events within 
this aperture. Scaling the background (495 counts in 
an annulus of 14,877 pixel area) to the source aperture, 
we found the average background contribution of 0.94 
counts. The background-subtracted, aperture-corrected 
source flux is Fx — (4.4±0.9) x 10~ 15 ergs s -1 cm~ 2 (the 
errors here and below are at the 68% confidence level for 
one interesting parameter), in the 0.4-2.0 keV band that 
includes the energies of all the 25 events detected (see 
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Fig. 2. — Energies and arrival times for 25 photons from the l'.'S 
radius aperture around B1257. 

Fig. 2). 

In Figure 2 we show the distribution of arrival times for 
the 25 detected photons over the duration of the obser- 
vation. The distribution of arrival times does not show 
any statistically significant deviations from the Poisson 
statistics. 

We fitted the spectrum in the 0.3-8 keV range with 
the absorbed BB and PL models. Since grouping of the 
detected 25 counts in energy bins would result in small 
numbers of counts per bin, using the standard \ 2 mini- 
mization technique would introduce a significant bias in 
the deduced model parameters (Cash 1979; Nousek & 
Shue 1989). Therefore, we have to use the C-statistic 
(implemented in XSPEC, ver. 11.3.2), without energy 
binning. To obtain constrained fits, we had to freeze 
the hydrogen column density, Ah- The pulsar's dis- 
persion measure, DM = 10.2 cm~ 3 pc, corresponds to 
JV H ~ 3 x 10 20 cm" 2 (assuming a 10% ISM ionization), 
which we adopt in our fits (see Table 1 and Fig. 3). 

The PL fit gives a photon index T = 2.75 ± 0.35 
and normalization constant Af — (1.7 ± 0.3) x 10~ 6 
photons cm~ 2 s _1 at 1 keV, corresponding to the lu- 
minosity L x = 4ird 2 F x nahs ~ (2.5 ± 0.5) x 10 29 d 2 00 
ergs s _1 in the 0.3-8 keV band, for isotropic emis- 
sion. The apparent temperature, T — 2.5 ± 0.3 MK, 
and the apparent projected area of the emitting region, 
Axa = 2.llo.g x 10 3 Gf| 00 m 2 , obtained from the BB fit 
are strongly correlated (see Fig. 3), which explains the 
large uncertainties of these parameters. Since Aj_ a is 
5 orders of magnitude smaller than the assumed pro- 
jected area of neutron star (NS) surface, ttR^ s , with 
-Rns ~ 10 4 rn 2 , such thermal radiation could originate 
only from small heated spots. If there are two identi- 
cal spots (polar caps) at star's magnetic poles, the ra- 
dius and the bolometric luminosity of each of the spots 
are R = (A Xa /fn) 1/2 ~ 26d 5 oo/~ 1/2 m and L hol = 
a(T a /g r )HA ±a /f) « 4.6 x lO 28 ^/- 1 ^ 4 ergs s" 1 , 
where the geometrical correction factor f < 1 depends 
on the angles ( (between the line of sight and the spin 
axis) and a (between the magnetic and spin axes) as well 
as the gravitational redshift factor g r = (1 — R g / i?Ns) 1 ^ 2 
(R g = 2953Mns/^o rn is the gravitational radius). Be- 
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Fig. 3.— Confidence contours (68% and 90%) for the PL (top) 
and BB (bottom) fits to the spectrum of B1257, for fixed JVjj val- 
ues (Nu = 0.3 X 10 21 cm -2 corresponds to the pulsar's dispersion 
measure, while Nn = 3 X 10 21 and 5 X 10 21 cm -2 were chosen to 
obtain X-ray efficiency of Bf257 similar to those of other pulsars; 
see §4.2). The PL normalization is in units of 10 -6 photons cm~ 2 
s — 1 keV - 1 at 1 keV. The BB normalization is the projected emit- 
ting area in units of m 2 , for d = 500 pc. The dashed lines are the 
lines of constant unabsorbed flux in the 0.3-8 keV band in the top 
panel, and constant bolometric flux in the bottom panel, both in 
units of 10 -15 ergs cm" 2 s _1 . 

cause of a substantial scatter and a shallow swing of the 
polarization position angle within the pulse profile in ra- 
dio polarimetric observations (Xilouris et al. 1998), the 
values of ( and a remain highly uncertain for B1257. 
However, thanks to the GR effect of bending photon 
trajectories in the NS gravitational field, the factor / 
varies in a relatively narrow range: 2R g /R < / < 1 (or 
0-83 < / < 1 for M NS = IAMq and R NS = 10 4 m, 
i.e. g r = 0.766), where the lower limit (corresponding to 
a = 0, C = 90° or C = 0, a = 90°) is estimated in the 
approximation outlined in Appendix of Zavlin, Shibanov, 
& Pavlov (1995). The BB model fits the spectrum better 
than the PL model (the C-statistic values are 77 and 84, 
respectively, for 524 degrees of freedom), but the PL fit 
cannot be rejected based on statistical arguments. 



4. DISCUSSION 

X-ray emission from a solitary MSP (or from an MSP 
in a wide binary) can be produced in the pulsar's magnc- 
tosphcre and/or hot polar caps, while an additional non- 
thermal component, emission from an intrabinary shock, 
can become important if the MSP is in a close binary sys- 
tem. The thermal (polar cap) and nonthermal (magneto- 
spheric and shock) contributions can be distinguished by 
the shape of the X-ray spectrum. If the observed emis- 
sion is predominantly nonthermal, the spectrum is best 
described by a PL, which is particularly hard when mag- 
netospheric emission dominates (e.g., PSR J0218+4232, 
B1821-24, B1937+20, for which L = 1-2; see Table 2). 
On the contrary, the spectrum of the thermal polar cap 
emission is quite soft; when fit with a PL model, it shows 
substantially larger photon indices, T = 3-5. Thermal 
emission is not seen in younger (r < 0.5 Gyr), power- 
ful (E > 10 35 ergs s _1 ) MSPs, while it dominates in 
solitary MSPs with E < 10 35 ergs s _1 . Observations of 
nearby thermally-emitting MSPs, such as PSR J0437- 
4715 (Zavlin et al. 2002; Zavlin 2006a,b), actually show 
both thermal and nonthermal components, the former 
dominating at lower energies, < 2 keV. 

Since the spin-down power of B1257 is < 10 34 erg s _1 , 
and its spectrum is quite soft (r ~ 3) and detected only 
below 2 keV, one can expect that the detected emission 
is mostly thermal emission from polar caps. If this is 
the case, the best-fit BB temperature, kT a « 0.2 keV, 
is about the same as found in other MSPs (e.g., PSR 
J0437-4715, J2124-3358, and J1024-0719). Although 
the projected emitting 2 x 10 3 m 2 , is rather 

uncertain, it is much smaller than the conventional po- 
lar cap area, A pc = 2Tr 2 R^ s /(Pc) = 1.05 x 10 7 m 2 
(for i?NS = 10 4 m) predicted by simple pulsar models, 
and, correspondingly, the best-fit BB polar cap radius, 
R ~ 30 m, is much smaller than the predicted radius, 
R pc = 1800 m ~ 60i?. Similar (albeit not so strong) dis- 
crepancies have been found for BB fits of other pulsar 
spectra, including ordinary old pulsars possibly showing 
thermal polar cap emission (e.g., Zavlin & Pavlov 2004; 
Kargaltsev, Pavlov, & Garmire 2006a). For many pul- 
sars, the discrepancy can be alleviated assuming that the 
polar cap is covered by a hydrogen or helium atmosphere, 
in which case the effective temperature would be a factor 
of 2 lower, and the radius a factor of 3-10 larger (Zavlin, 
Pavlov, & Shibanov 1996). However, the discrepancy is 
too high for B1257 to be explained in this way. An- 
other explanation for such a discrepancy was suggested 
by Zhang, Sanwal, & Pavlov (2005), who analyzed the X- 
ray emission from the ordinary drifting pulsar B0943+16 
and suggested that only a small fraction of the polar cap 
surface, associated with footprints of sparks produced by 
intermittent breakdowns of an "inner gap" above the po- 
lar cap, is hot enough to emit X-rays. However, the "fill- 
ing factor", A/ A pc ~ 2 x 10~ 4 , is a factor of 100 smaller 
for B1257 than for B0943+16, and it is currently unclear 
whether or not this hypothesis is applicable to MSPs. 
The area would become much larger if the actual Ah 
is substantially larger than the 3 x 10 20 cm -2 estimated 
from the pulsar's dispersion measure (e.g., A±, a <~ 2x 10 5 
m 2 for Ah = 3x 10 21 cm~ 2 , as demonstrated in Fig. 3), 
but such large Ah values would require an additional ab- 
sorber, such as a dust cloud near the pulsar (see §4.2), 
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because they strongly exceed the total Galactic TVh in the 
direction to B1257 (e.g., N m = 2.2 x 10 20 cuT 2 ; Dickey 
& Lockman 1990). 

The magnetospheric interpretation of the X-ray emis- 
sion from B1257 looks less plausible than the thermal 
one, but we cannot firmly rule it out with the small num- 
ber of counts detected. Although the best-fit power-law 
is unusually soft for magnetospheric emission, the uncer- 
tainty of photon index allows smaller values of T (see 
Fig. 3), possibly compatible with the magnetospheric in- 
terpretation. 

4.1. Intrinsically underluminous millisecond pulsar? 

Whether the detected X-ray emission is thermal or 
nonthermal, the B1257's apparent luminosity is very 
low in comparison with other MSPs detected in X- 
rays, and its apparent X-ray efficiency 2 , rj = Lx/ E « 
1.3 x 10~ 5 d 2 00 (l - 0.59d 50 o)~\ is smaller than the typ- 
ical values, i] ~ 10~ 4 -10~ 2 - 5 . This is demonstrated in 
Figure 4, which shows the 0.1-10 keV luminosities and 
spin-down powers for 31 MSPs (we chose the 0.1-10 keV 
band following Bogdanov ct al. 2006, whose results on 
14 MSPs in the globular cluster 47 Tuc are included in 
Fig. 4). The X-ray luminosities and the corrections of 
spin-down power for the Shklovskii effect (see §1) were 
calculated using the distances and transverse velocities 
listed in Table 2. To make the picture more legible, we 
chose not to show the Lx and E uncertainties, except 
for B1257 (a comprehensive discussion on the Lx uncer- 
tainties is presented by Bogdanov et al. 2006). For the 
MSPs showing predominantly nonthermal X-ray emis- 
sion (blue squares and green triangles in Fig. 4), we plot 
the "isotropic luminosities", Lx = 4Trd 2 F x lDahs , which 
can be higher or lower than the actual luminosities, de- 
pending on the (unknown) intrinsic angular distribution 
of the pulsar's radiation. For the MSPs whose X-ray 
emission is most likely thermal (red circles in Fig. 4), 
we used the "equivalent sphere luminosities" (corrected 
for the general relativity effects assuming g r = 0.766) 
instead of the "true" polar cap luminosities, because 
the angles £ and a needed for calculation of the lat- 
tcrs (see §3) are not known for most of the MSPs. For 
instance, the equivalent sphere bolomctric BB luminos- 
ity is L b s ol = A(A ±a g 2 r )a(T a /g r ) i = (2/ 5 2 )(2L bol ), where 
2Lboi is the luminosity of two polar caps. The factor 2/<? 2 
is not very different from 1 for any £ and a at the ex- 
pected values of g r ; e.g., it varies between 0.97 (at ( = 0, 
a = 90°) and 1.17 (at C = a = 0) at g r = 0.766. 

Since both Lx and E depend on distance, which has 
not been accurately measured for B1257, the measured 
X-ray flux (with account for measurement uncertainties) 
corresponds to a strip in the E-Lx plane, shown in Fig- 
ure 4 for the d = 300-800 pc range (the dashed curve 
corresponds to the best fit). We see that for the most 
plausible distances, d < 750 pc, not only the B1257's 
X-ray luminosity is lower than for any MSP detected 
in X-rays, but also its X-ray efficiency is lower, except 
may be PSR J0034-0534 (#2 in Fig.' 4), for which we 

2 We emphasize that this low efficiency is obtained assuming 
Nn = 3 X 10 20 cm~ 2 , estimated from the pulsar's dispersion mea- 
sure. The inferred luminosity and efficiency would be higher if a 
larger TVh value is assumed, but they still would be lower than 
typical ones as long as the X-ray absorption is caused by the ISM. 



know only an upper limit on E because its proper motion 
has not been measured. Only for distances approaching 
the limiting value of 848 pc (see §1), rj becomes similar 
to those of the majority of MSPs (but the luminosity, 
Lx —* 9 x 10 29 ergs s _1 , still remains lower than those 
of most of the other MSPs). If the observed B1257's ra- 
diation is magnetospheric, the low apparent luminosity 
might be explained by an unfavorable direction of the 
pulsar beams (i.e., the actual luminosity of B1257 being 
larger than the isotropic luminosity). However, it cannot 
be explained under the more plausible assumption that 
it is thermal radiation from two symmetric, isotropically 
emitting polar caps. Under this assumption, the maxi- 
mal value of 2Lboi/^boi = (2.fffr) _1 does not exceed 1.03 
for g r — 0.766 (see above), i.e. the maximum luminos- 
ity of two polar caps corresponding to the observed flux 
is only 3% higher than the equivalent sphere luminosity 
plotted in Figure 4. To get a larger intrinsic luminosity, 
one might speculate that the caps are very asymmetric 
(e.g., one of them is much brighter than the other and 
is invisible for most part of pulsar period, which can oc- 
cur only when both Q and a are very small. Also, the 
polar cap emission can be anisotropic, contrary to BB 
emission, because of the limb-darkening effect in the NS 
atmosphere, but this anisotropy is not strong at the rel- 
atively low magnetic fields of MSPs (see Zavlin et al. 
1996). Thus, although there are ways to derive larger 
luminosity and efficiency from the same observed flux, 
these quantities arc not expected to be much larger than 
our current estimates, at least if the observed emission is 
indeed thermal. 

Among the MSPs detected in X-rays, of particular in- 
terest is PSR B1620-26 (#11 in Fig. 4), a member of a 
triple system that likely contains a Jupiter-mass planet 
on a distant orbit around the inner pulsar + white dwarf 
binary (F b in = 191 d, P p i a not ~ 100 yr; Sigurdsson & 
Thorsett 2005, and references therein). The globular 
cluster M4, which hosts PSR B1620-26, was observed 
with Chandra ACIS on 2000 June 30 for 25.8 ks (Bassa et 
al. 2004). Since the description of the results on B1620- 
26 is very sketchy in that paper, we reanalyzed the data 
using the same approach as for B1257 (§3). The pul- 
sar was detected 1'.2 off- axis, and its image looks slightly 
extended in the north-south direction. In an elliptical 
aperture with major and minor axes of 2'.' 7 and l'/7 we 
found 21 photons in the 0.3-8 keV band and measured 
the flux F x = (4.2 ± 1.0) x 10~ 15 ergs cm~ 2 s" 1 . We fit 
the source spectrum with the PL and BB models at fixed 
iV H = 2.36 x 10 21 cm" 2 . The PL fit gives r = 2.3 ± 0.4 
and Fx nabs « 7.7 x 10~ 15 ergs cm~ 2 s^ 1 , correspond- 
ing to the isotropic luminosity of 2.8 x 10 30 ergs s _1 in 
the 0.3-8 keV band (L x = 4.6 x 10 30 ergs s" 1 in the 
0.1-10 keV band) at the distance of 1.73 kpc to the glob- 
ular cluster. The BB fit gives kT a = 0.45 ± 0.08 keV, 
A ±3 , = l-0±g;| x 10 3 m 2 , and L b s ol « 3.0 x 10 30 ergs 
s^ 1 . Similar to B1257, both fits are statistically accept- 
able and give about the same X-ray luminosity. The 
apparent BB temperature is higher than for the other 
thermally emitting MSPs, which indicates that we possi- 
bly detect both thermal and nonthermal radiation, per- 
haps even some contribution from a pulsar wind neb- 
ula, as indicated by the possible extension of the source 
image. Regardless of the emission mechanism, the lu- 
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minosty of this MSP is substantially higher than that 
of B1257. Unfortunately, we do not know its intrin- 
sic spin-down power (hence, X-ray efficiency) because of 
poor knowledge of parameters of the putative planet that 
strongly contributes to the observed value of P. In Fig- 
ure 4, we plotted the point corresponding to B1620-26 at 
E = 0.01-Eobs, following the assumption by Thorsett et 
al. (1999), and showed the E uncertainty by the double- 
sided horizontal arrow. Notice that if E is indeed so 
small, then B1620-26 is the most efficient X-ray emit- 
ter (rj ~ 10~ 2 ) among the 31 MSPs, opposite to the 
other planet pulsar B1257. We should remember, how- 
ever, that the planet in the B1620— 26 system is quite 
different in properties and history from the B1257 plan- 
etary system. 

4.2. Absorbed by matter orbiting the pulsar? 

The planets around B1257 were likely formed from a 
disrupted or ablated stellar companion that had possi- 
bly provided the matter to spin up the pulsar to its mil- 
lisecond period (Phinney & Hansen 1993; Podsiadlowski 
1993). Some material left over the planet formation (as- 
teroids, meteoroids, dust) can still rotate around the pul- 
sar (e.g., Bryden et al. 2006; Cordes & Shannon 2006). 
Absorption of the B1257's radiation by such circumpusar 
material might explain the relatively low observed X-ray 
flux. 

To explore this possibility, we assume that X-rays pass 
through a cloud that contains some grains or rocks. For 
optically thick grains (radius a lpm at E <~ 1 keV), 
the grain cross section, a g <~ na 2 , is independent of pho- 
ton energy, so the absorption by grains does not affect the 
shape of X-ray spectrum. To obtain the intrinsic B1257's 
luminosity and efficiency similar to those of other MSPs, 
the optical depth r of the intervening cloud should be in 
a range of 2-5, along the line of sight, corresponding to 
the column density N g ~ r(7ra 2 ) -1 and number density 
n 9 ~2x 10 _14 ra _2 Z _1 cm~ 3 , where a is the grain radius 
in units of cm, and I is is the path length through the 
cloud in AU. The mass of the putative cloud depends 
on its geometry and location with respect to the line of 
sight, as well as on size and composition of grains, all of 
which are unknown. It can be scaled as 

M cl ~ 3 x 10 26 T(a/l)pV cl g , (1) 

where p is the mass density of the grain material in 
gem -3 , and V^i the cloud's volume in units of AU 3 . 
Thus, at t ~ 3, p ~ 2 g cm -3 , and a characteristic 
cloud size of 1 AU, we obtain M c] - (0.003-300)M e 
for a — 0.01-1000 cm. If the absorption of X-rays is 
caused by large bodies, the mass of the cloud becomes 
uncomfortably large (e.g., in comparison with the total 
mass of B1257's planets, w 8M®), especially if the mass 
distribution is concentrated towards the planets orbital 
plane. Since the orbital plane is inclined by 40° to the 
line of sight (Konacki & Wolszczan 2003), only a periph- 
ery of such a distribution (a flaring disk) can contribute 
to the absorption, and this would imply V c \ ^> 1. On 
the other hand, distant orbits of circumpulsar matter 
can be strongly inclined with respect to the planets or- 
bital plane, so that one might imagine an azimuthally 
nonuniform belt of matter on an orbit partially eclips- 
ing the pulsar. An argument against this hypothesis 



is that grains/rocks are expected to have been evapo- 
rated/ablated by the pulsar wind unless they are very 
large (e.g., > 100 m in size; see Miller & Hamilton 2001) 
or replenished by collisions of larger bodies with a rate 
exceeding the rate of evaporation/ablation. 

Absorption of X-rays by very small, optically thin 
grains is virtually indistinguishable from absorption by 
the ISM gas (Wilms et al. 2000, and references therein). 
Therefore, its effect on the absorbed spectrum can be 
crudely modeled by increasing Ah in the ISM absorption 
models. To explore this possibility, we fitted the X-ray 
spectrum of B1257 with the same PL and BB models 
as previously but assuming larger (fixed) values for Ah- 
We show the parameter confidence contours in Figure 3 
for N H = 3 and 5 x 10 21 cur 2 for the PL and BB fits 
respectively (these values provide a factor of ~10 higher 
luminosities for the two models) . The PL fit yields a very 
large photon index, T ~ 4-5, suggesting that the spec- 
trum is, in fact, thermal. The BB fit gives a slightly lower 
temperature, T a = 1.51 ±0.16 MK and a much larger 
emitting area, A ±a = 1.9±gJ x 10 5 m 2 (~ 0.02A pc ), cor- 
responding to a luminosity Lboi ~ 5.6 x 10 d§ o/ ~ 1 g^ i 
ergs s _1 . A factor of 10 higher X-ray efficiency obtained 
from this fit is similar to those of most of MSPs. The 
larger Ah required for this increase can be used to esti- 
mate the mass of the intervening cloud, 

M cl ^2xl0 24 eU c i/- 1 g (2) 

where £ ~ Z st /Z is a factor depending on element abun- 
dances in the grains (Z st « 0.02 is the 'standard' metal- 
licity used in the absorption model, and Z is the actual 
metallicity in the grains). We see that a much lower 
mass, M c \ ~ 3 x 1O~ 4 £M0 for a 1 AU characteristic 
cloud size, is needed to explain the low X-ray efficiency of 
B1257 by absorption in a circumpulsar dust of micropar- 
ticles. However, such small particles could be blown out 
from the pulsar's vicinity by the pulsar wind and radia- 
tive pressure (especially at earlier epochs when the spin- 
down liminosity was higher), and they could be swept 
out from the outskirts of the planetary system by the 
ram pressure of the oncomimg medium. Therefore, to 
explain the low apparent X-ray luminosity and efficiency 
as due to absorption by small particles, we have to as- 
sume that the particles are being replenished by collisions 
of larger bodies in the putative asteroid belt. 

Some additional information on the circumpulsar mat- 
ter can be provided by infrared (IR) observations. Since 
such matter is being heated by the pulsar's radiation 
(photons and pulsar wind), it should emit IR radia- 
tion whose spectrum depends on grain's temperature and 
composition. So far, searches for such emission at A ~ 
10-1000 pm have yielded only upper limits (Bryden et 
al. 2006, and references therein), which, however, can be 
used to put some constraints on the circumpulsar matter 
properties. For instance, from the upper limit on spec- 
tral flux, F v < 45 pJy at A = 24 pm (Bryden et al. 2006), 
we obtain an upper limit, A < 1.5 x 10 23 [exp(600 K/T) - 
1]?24 ^ioo cm2 j on the emitting area of grains in the cloud 
(^24 is the emission efficiency q\ at A — 24 /im). If the 
cloud is optically thin in IR, this limit translates into 
M cl - Aap/3 < 0.5 x 10 23 ap [exp(600 K/T) - l]^ 1 ^ 
g (e.g., M c i < 4 x 10 25 a g for T = 100 K, p = 2 g 
cm~ 3 , a > 0.5 x 10~ 2 cm), and, together with the es- 
timates of cloud mass needed to explain the additional 
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absorption (equations [1] or [2]), it constrains the cloud 
size. For instance, in the case of large grains, we ob- 
tain V c i/l < 1.7 x l(T 4 [cxp(600/T) - l]^ 1 AU 2 (e.g., 
V c i/l < 2.2 x 1CT 2 AU 2 for T = 100 K, r = 3). Since 
such estimates strongly depend on the unknown temper- 
ature, it would be very useful to obtain better constraints 
on the temperature and emitting area from deeper IR ob- 
servations. 

Overall, the hypothesis that the low apparent X-ray ef- 
ficiency is caused by absorption of X-rays in circumpulsar 
matter does not look unreasonable at this point. How- 
ever, it implies a rather narrow range of optical depths 
(and, in the case of large grains, rather large masses of 
absorbing matter) to make the B1257's X-ray efficiency 



consistent with those of other MSPs. In addition to deep 
IR observations, a possible way to confirm or reject this 
interpretation would be monitoring of the pulsar's X-ray 
emission. If the X-ray flux shows substantial variations 
(e.g., on a year timescale), it could be caused by vari- 
able absorption in an orbiting, nonuniformly distributed 
circumpulsar matter. 
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Fig. 4. — X-ray luminosty in the 0.1-10 keV band versus intrinsic spin-down powe r for 31 MSPs observed with Chandra and XMM- 
Newton. The numerical labels c orrespond to the numbers listed in first column of Table 2, while the letters mark the MSPs in 47 Tuc 
(Table 4 in Bogdanov et al. 2006). The red circles, blue squares, and green triangles correspond to MSPs whose X-ray emission is dom 
inated by thermal polar cap, magnetospheric, and intrabinary shock emission, respectively. The open and filled symbols denote solitary 
and binary pulsars, respectively. MSPs in globular clusters are marked by black + signs. The double-sided horizontal arrow attached to 
#11 (PSR B1620-26 in M4) demonstrates the uncertainty of its intrinsic E (see text). For #2 (PSR J0034-0534), only an upper limit on 
E is known because the lack of proper motion measurements makes the correction for the Shklovskii effect impossible. The curved strip 
shows possible positions for B1257 in the E-L^ plane for a range of distances, 300—800 pc (the dashed curve within the strip corresponds 
to the best-fit luminosity for a given distance, and the point #9 corresponds to d = 500 pc). Three dashed straight lines correspond to 
constant values of the X-ray efficiency, n = Lx/E. 
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TABLE 1 

Fitting parameters for the PL and BB models 



Model 


Nh,20 


A/" a or A ±a b 


T or kT c 


C/dof 


L x or ^bol d 


PL 
BB 
PL 
BB 


3 
3 

30 
50 


i 79 +U.34 
1 - ' z -0.34 

2 1+ 1 ' 9 
4 QO+0.84 

lqn +260 
iyu -l()() 


z - '°-0.36 
n 91C+O.O25 
U.z±o_ 023 

4 47+0.46 

o no+o oi5 

U.ldU_ ()14 


84/524 
77/524 
77/524 
77/524 


2 4 7 +U.5U 

z -^'-0.48 
J-O4-0.40 

16.615;? 

99 q+9.2 



Note. — The fits arc for fixed Nu.2t) = JVh/10 20 cm 2 (second column). 
a Spectral flux in units of 10 -6 photons cm -2 s _1 kcV -1 at 1 keV. 
b Projected area of emitting region for the BB model (in 10 3 m 2 ) for d = 500 pc 
C BB temperature in keV. 

d Unabsorbed luminosity in the 0.3-8 keV band or apparent bolometric luminosity (iaT^A± a ), in units of 10 29 ergs s _1 for d = 500 pc. 



TABLE 2 

Millisecond pulsars observed with Chandra and XMM-Newton 



# 


Name a 


Type b 


P 


P C 

-*orb 


v ± d 


log E c 


S 


rs 


logL x h 


- log rf 


Refs. 








ms 


d 


km/s 




kpc 










1 


J0030+0451 


th 


4.87 




<16 


33.53 


0.30 J 


4.7 


30.60 


2.93 


1,2 


2 


J0034-0534 k 


th? 


1.88 


1.59 




<34.48 


0.53 


2.5 


30.00 


<4.48 


3,4 


3 


J0218+4232 


mag 


2.32 


2.03 


50 


35.38 


2.6 


1.1 


32.53 


2.85 


5,6 


4 


J0437-4715 


th 


5.76 


5.74 


106 


33.46 


0.16 j 


4.1 


30.72 


2.74 


7,8 


5 


J0737-3039A 1 


mag? 


22.70 


0.10 


10 


33.76 


0.48 


3.2 


30.30 


3.64 


9-12 


6 


J0751+1807 


sh 


3.48 


0.26 


32 


33.86 


1.12 


1.6 


30.83 


3.03 


13,14 


7 


J1012+5307 


sh? 


5.26 


0.60 


102 


33.43 


0.84 m 


1.8 


31.04 


2.39 


13-16 


8 


J1024-0719 


th 


5.16 




109 


32.63 


0.39 


3.7 


29.88 


2.75 


7,8 


9 


B1257+12 


th 


6.22 




224 


33.89 


~0.50 


2.7 


29.50 


4.44 


17,18 


10 


B1534+12 1 


sh? 


37.90 


0.42 


120 


33.23 


l.Qi 


3.2 


30.15 


3.08 


10,19 


11 


B1620-26 n (M4) 


th? 


11.07 


191 




<34.20 


1.73P 


2.4 


30.48 


<3.72 


20,18 


12 


J1740-5340 (NGC 6397) 


sh 


3.60 


1.34 




34.52 


2.55 p 


1.5 


31.30 


3.22 


21,22 


13 


B1821-24 (M28) 


mag 


3.05 




120 


36.35 


5.5? 


1.2 


33.22 


3.13 


23 


14 


J1911-6000C (NGC 6752) 


th? 


5.28 






32.77 


4. 1? 


2.5 


30.64 


2.13 


24 


15 


B1937+21 


mag 


1.56 




14 


36.04 


3.55 


1.9 


33.18 


2.86 


25 


16 


B1957+20 


sh 


1.61 


0.38 


359 


34.87 


2.49 


1.9 


31.82 


3.05 


26,27 


17 


J2124-3358 


th 


4.93 




62 


33.63 


0.27 


3.3 


30.26 


3.37 


4,28 



Note. — References: 1- Becker & Aschcnbach 2002, 2- Lommen ct al. 2006, 3- Bailcs et al. 1994, 4- Zavlin 2006a, 5- Navarro et al. 
1995, 6-Webb et al. 2004a, 7- Hotan et al. 2006, 8- Zavlin et al. 2002, 9- Lyne, et al. 2004, 10- Kargaltsev et al. 2006b, 11- Kramer et al. 
2006, 12- Chatterjee et al. 2007, 13-Webb et al. 2004b, 14- Nice et al. 2005, 15- Lange et al. 2001, 16- Callanan et al. 1998, 17- Wolszczan 
et al. 2000, 18- this work, 19- Stairs et al. 2002, 20- Bassa et al. 2004, 21- Bassa & Stappers 2004, 22-Grindlay et al. 2002, 23- Becker et 
al. 2003, 24- D'Amico et al. 2002, 25- Nicastro et al. 2004, 26- Toscano ct al. 1999b, 27- Stappers et al. 2003, 28- Toscano et al. 1999a. 
a For the globular cluster MSPs, the host cluster is given in parentheses. 

b Type of emission dominating in the X-ray range: thermal (th), magnctospheric (mag), or emission from an unresolved shock (sh). The 
most uncertain cases are marked with '?'. 
c Orbital period for binary MSPs. 
d Transverse velocity. 

intrinsic spin-down power, corrected for the Shklovskii effect and the effect of gravitational potential for globular cluster pulsars. 
'Distances estimated from the model for Galactic electron distribution by Cordes & Lazio (2002) unless indicated otherwise. 
g Photon index from fitting with one-component PL model (irrespective of fit quality), characterizing the spectral hardness. 
h Best estimate for the unabsorbed X-ray luminosity in the 0.1-10 keV band. 
l n = Lx/E is the total X-ray efficiency. 
J Distance measured from radio-timing parallax. 

k Since proper motion was not measured for this pulsar, we cannot correct E for the Shklovskii effect. 
'Double neutron star binary. 
m Distancc estimated from observations of the white dwarf companion. 

"Triple system with a planet. The intrinsic P and E are constrained very poorly (see text). 
p Distancc to the globular cluster. 



